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For  an  isotropic  solid  shell,  a yield  condition,  a strain-hardening  rule, 
and  a flow  rule  have  been  formulated  in  terms  of  the  membrane  forces  and 
bending  moments  and  the  strains  and  curvatures  of  the  middle  surface.  The 
results  of  the  proposed  theory  are  '•ompared  to  the  results  of  the  through- 
the-thickness  integration  of  the  stress-strain  relations  of  the  material. 
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LIST  OF  SYMBOLS 


parameters  in  the  equation  of  the  .imit 
yield  surface  Eq . (13) 

= shell  thickness  (solid  shell,  Fig  2) 

= dimensions  of  a sandwich  shell  (Fig.  2) 

= indices,  i = 1,  2;  j = 1,  2 
= strains  components  of  the  middle  surface 
= curvature  components  of  the  middle  surface 

= see  Fig.  2 

= column  matrix  of  shell  strains  (Eq.  (22)) 

= plastic  part  of  e 

= column  matrix  of  stress  resi'Itan  ts , (Eq . (21)) 

= time 

= expression  defined  by  Eq . (31) 

= elasto-plasf.ic  tangent  stiffness  matrix 
= elastic  moduli  matrix 

" subsequent,  initial,  and  limit  yield  functions 

- absolute  values  of  the  gradients  of  the  yield 
function  F (see  Eqs.  (IS)  and  (19)) 

= resultant  stress  invariants  (see  Eqs.  (6),  (7),  (8)) 

membrane  forces  in  shell 
moments  in  shell 
strain  hardening  parameters 

for  solid  shells 

for  ideal  sandwich  shells 
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LIST  OF  SYMBOLS  (Continued) 

a = a parameter  in  the  yield  condition 

3 = a parameter  in  the  hardening  rule 

X = a parameter  in  the  flow  rule 

v = Poisson  ratio 

a..  = stress  components 

iJ 

Oq  = yield  stress  in  uniaxial  tension 
£ . . = strain  components 
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INTRODUCTION 


This  report  discusses  the  elasto-plastic  behavior  of 
thin  plates  and  shells  whose  material  is  an  elasto-plastic 
solid,  linear  in  the  elastic  range,  obeying  Mises'  yield 
condition,  and  deforming  plastically  according  to  the 
Prand tl-Reuss  flow  rule.  The  mechanical  properties  of  this 
solid  are  characterized  by  two  elastic  constants  and  the 
yield  stress  in  uniaxial  tension,  0^. 

The  plate  or  shell  is  assumed  to  be  a solid  layer,  with 
thickness  h;  occasionally,  reference  will  be  made  to  an  ideal 
sandwich  shell  (Fig.  1)  such  that  d >>  t.  The  basis  assumption 
of  the  plate  and  shell  theory  used  in  this  paper  are  summarized 
in  the  expressions  for  the  strain  components,  at  any  point  of 
the  shell  in  terms  of  the  strain  and  curvature  components  of 
the  middle  surface 


(z) 


e . . + k 
id 


(1) 


and  in  the  expressions  for  the  membrane  forces  and  the  bending 

moments  (stress  resultants) 

fh/2  [h/2 
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-h/2 


0 . . d: 
13 


M, 


13 


a . . z d z 
13 

-h/2 


(2) 


It  becomes  evident  from  the  above  equations  that,  within  the 
accepted  approximation,  there  is  no  distinction  between  a 
plate  and  a shell,  as  far  as  the  stress-strain  relations  are 
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concerned,  and  the  term  "shell"  will  be  used  to  describe 
the  structure  under  consideration. 


The  strain  - normal  force  and  the  curvature  - moment 
relations  under  the  conditions  of  uniaxial  stress  (i.e.  beam 
behavior)  are  illustrated  in  Fig.  2.  The  specific  objective 
of  this  work  is  the  development  of  the  relations  between  the 
stress  resultants,  N„  and  , and  the  strain  components, 
e_  and  k^  , for  a general  type  of  loading. 

The  yield  conditions  for  the  simplified  behavior  in 
bending  have  been  proposed  in  several  previous  investigations; 
a critical  evaluation  and  comparison  of  the  existing  yield 
surfaces  can  be  found  in  a paper  by  M.  Rob inson  (Ref . [1]).  In 

order  to  take  into  account  the  actual  moment  - curvature 
relation,  M . A . Crisfield  (Ref.  [2]),  applies  the  ideas  of 
isotropic  strain  hardening  of  the  theory  of  plasticity  to  the 
elas to-plas tic  behavior  of  shells.  The  approach  reported  in 
this  paper  is  a continuation  and  extension  of  the  above 
earlier  works. 


The  fact  that  the  stress  resultants  N..  and  MJ . of  the 

ij  ij 

classical  shell  theory  are  not  sufficient  to  describe  the 
state  of  stress  has  been  recognized  by  G.  Wempner.  In  Ref.  [3], 
he  introduced  certain  higher  - order  moments  which,  together 
with  the  classical  stress  resultants,  form  the  dynamic  variables 
of  the  problem.  In  Wempner's  most  recert  work,  described  in 
a private  communication  to  this  author,  the  stress  distribution 
through  the  thickness  of  the  shell  is  expanded  in  teims  of 


( 
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Legendre  polynomials,  with  the  coefficients  of  this  expansion 
being  the  additional  dynamical  variables. 

It  is,  of  course,  possible  (and  it  has  been  dene  in 
many  investigations  of  elasto-plastic  shells)  to  avoid 
completely  the  problem  of  the  shell  constitutive  equations. 

In  the  "through-the-thickness-integration"  approach,  for 
given  increments  of  e..  and  k..,  the  increments  of  strains 
e.j(z)  are  determined  with  Eqs.  (1);  then,  from  appropriate 
constitutive  equations  of  the  material,  the  increments  of 
stresses  0^(2)  are  computed;  finally,  the  increments  of  shell 
forces  and  moment,  are  determined  by  numerical  evaluation  of 
the  integrals  in  Eqs.  (2).  Although  workable  and  accurate, 
the  procedure  requires  sometimes  prohibitively  large  storage 
capacity  of  the  computer. 
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II  INITIAL  AND  LIMIT  YIELD  SURFACES 


The  initial  yield  condition,  or  the  initial  yield  surface 
in  the  stress  space,  can  be  easily  derived  from  our  basic 
assumptions.  The  stress  components  at  the  top  and  bottom 
surfaces  of  a solid  shell  are 

N . . 6 M . . 

a . . = -jp-  + — ip- 

13  h - .2 


(3) 


where  the  plus  sign  applies  to  the  top  and  the 

minus  sign  to  the  bottom  of  the  shell.  Substitution  of 

expression  (3)  into  Mises'  yield  condition 


2taU  + °22  “ °11  °22  + 3 °12)  * 1 

0 


results  in 


where 


J -hi  + 21  =1 

N M - NM 


12  2 2 

I = — — (N  +N  -N  N +3N  ) 

N 2'  11  22  N11  22  12; 

K0 

12  2 2 

I = — =—  (M  + M - M M + 3 M ) 

2k 11  22  11  22  12' 


M 


0 


(A) 

(5) 

(6) 

(7) 


NM 


N0M0v‘'ll  “11  ' “22“22  ” 2 N11M22  ~ 2 N22M11  + 3 N12M12^  ^8) 


*(N  M + N M 


4.'  ’ 


No  = aoh 


Mo  = V /6 


(9) 


When  the  expression  (5)  is  used,  the  top  or  the  bottom  of 
the  shell  should  be  considered  in  order  to  have  the  larger, 

[NM  ' 


i.e.,  positive,  value  for  + 2 I . This  is  assured 
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by  writing  the  initial  yield  condition  of  a solid  shell  as 


F0  5 *»  + I*  + 2 Urn.  I - 1 


N 


M 


NM1 


(10) 


In  tv'  case  of  an  ideal  sandwich  shell  the  stresses  are 
computed  from 


id 


N . . M . . 
_i2  + _il 
2 1 - dt 


and  an  identical  argument  leads  to  the  condition  (10)  except 
that  Nq  and  1.^  are  now 


0O  2 1 


M0  = 0Q  td 


The  physical  meaning  of  the  Quantities  and  MQ  is 
as  in  Fig.  2.  The  corresponding  strains  are  e^  and  k^, 


respectively.  For  a solid  shell 

0, 


0 


'0 


k0 


V 

Eh 


(11) 


(12) 


No  direct  derivation,  of  the  type  used  above,  is  possible 
for  the  limit  surface.  Instead,  approaches  which  are  essentially 
surface-fitting  procedures  have  proven  to  be  useful.  Suppose 
that  the  limit  surface  is  represented  by  a linear  equation  in 

IN’  IM’  and  1NM 


FL  = 3 IN  + b IM  + C XNM  ~ 1 


In  principle,  the  parameters  a,  b,  and  c should  be 
determined  in  such  a way  as  to  minimize  tha  difference  between 
the  surface  (13)  and  the  yield  surlace  determined  by  some  more 
precise  calculations.  For  practical  purposes,  however,  the 


(13) 
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following  argument  results  in  sufficiently  accurate  values 

of  a,  b,  and  c.  By  considering  the  case  of  membrane  forces 

only,  we  find  that  with  a = 1 Eq . (13)  becomes  the  exact 

2 2 

limit  condition.  Similarly  with  b = /M  , Eq . (13)  will 
produce  the  exact  results  for  the  cases  of  bending  moments 
only.  It  is  now  necessary  to  consider  one  more  loading  case, 
with  known  exact  solution,  in  order  to  find  the  value  of  c. 
Such  a case  may  be  taken  as  corresponding  to  the  maximum 
value  of  I . Upon  reflection,  it  becomes  evident  that  it 

iS  N11  = N22  » M11  = M22’  N12  = °’  and  M12  = °- 

The  stress  distribution  in  a section  is  then  as  shown 
in  Fig.  3,  and  the  maximum  of  corresponds  to  n = h/2/3. 

A simple  computation  results  in 

h,  ■ j • hi  - «i2/9ko2>  htf  - 2 ^ Ml/9M0 

Equation  (13)  will  be  satisfied  with  the  above  values  if 
c = Mq/Ml  /3  . 

This  form  of  F is  identical  with  Iliushin's,  (Ref.  [4]) 

JL 

yield  condition  obtained  from  different  considerations. 

For  an  ideal  sandwich  shell,  the  functions  F and  F 

U L 

coincide,  since  there  is  no  distinction  between  initial 
yield  and  limit  state. 


Ill  PROPOSED  YIELD  CONDITION  AND  HARDENING  RULE 


The  following  yield  condition  is  proposed  to  describe 
the  "subsequent"  yield  surfaces  as  the  loading  path  moves 
from  the  initial  yield  surface  towards  the  limit  surface 


F 5 + TM  + “ IV  ’ 1 


(14) 


where 


h = TT^ll  “ Mll)2  + <M22  - M22)2 

M0 


*s2 


^Mll  " M11^M22  " M22^  + 3^M12  “ M12^  ^ 


The  quantities  M^_.  , which  will  be  referred  to  as  "hardening 
parameters",  are  defined  by  the  following 


(15) 


SF  ' 9F 

If  F = 1 and  N . . + w— - M.  . > 0: 

3N,,  oM . . ij 


ij 


0 FM 


(16) 


If  F < 0 or  ——  N.  . + M.  . < 0: 

oN . . in  oM . . in  - 

in  id 


d M . . = 0 


(17) 


The  symbols  Fr  and  F^  are  defined  as 

2 2 


Fs  - [(N 


+ (M 


0 3n1]l^ 

* <N0 

) 

3N 

N22 

+ <«0 

3N 

3F  s.2 

•\\t  ) 

+ (M„  ■ 

3F  2 

*\  w ) 

+ (M„  ' 

3F 

12 

2 1/2 
-)  ] 


11 


22 


12 


(18) 


r 


14 


^ TT  2 ^r*  2 j\-ri  2 1/2 

“M0  35^  + (M0  3M^')  + <M0  3^  1 


F is  evidently  the  absolute  value  of  the  vector  grad  F,  in 
s 

a dimensionless  formulation;  F^  is  the  part  of  grad  F which 
corresponds  to  the  bending  moments  only. 


The  function  F which  appears  in  Eq . (16)  is  defined  by 

Li 

Eq . (13).  With  the  values  for  a,  b,  and  c as  determined  in 

the  preceding  section  and  with  M /M.  = 3/2  (for  solid  shells), 

L U 

it  reads 


FL  " XN  + 9 *M  + 


3 /T 


The  above  formulation  contains  two  parameters,  a and  3* 
The  parameter  a should  be  variable.  When  the  loading  path  is 
still  in  the  initial  elastic  range,  its  value  should  be 
a = 2 to  assure  correct  predictions  of  the  instant  of  first 
yielding.  As  the  loading  path  approaches  the  limit  surface, 
the  value  of  a should  approach  the  value  of  c in  Eq . (13). 

It  appears  however,  that  sufficiently  close  approximations 
to  the  exact  results  can  be  obtained  with  a constant  value 


of  a.  Here,  a = has  been  used,  i.e.  the  limit  surface 

3/3 

is  reproduced  correctly  while  an  error  is  accepted  in  the 
initial  yield  surface. 


The  parameter  3 controls  the  moment  - curvature  relation 
in  the  plastic  range.  Again,  a constant  value  of  3 = 2,  has 
been  found  reasonably  satisfactory  for  solid  shells. 
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The  hardening  law  represented  by  Eq . (16)  is  neither 


i 

i 


i 


v > 

1 


I 
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isotropic  nor  kinematic.  Its  choice  is  motivated  solely  by 
the  fact  that  it  reproduces  fairly  closely  the  actual 
behavior  of  a solid  shell  in  the  plastic  range.  It  reproduces 
also  the  lowered  yield  point  ("Bauschinger  effect")  which 
manifests  itself  if  the  bending  moment  is  reversed,  the  shell 
unloaded  and  then  loaded  in  opposite  direction. 
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IV  FLOW  RULE  AND  TANGENT  STIFFNESS 


To  complete  the  formulation  of  the  behavior  of  eiasto- 
plastic  shells,  it  is  necessary  to  state  the  elastic  law 
and  the  flow  rule.  For  this  purpose,  the  stress  resultants 
and  the  strain  components  of  the  shell  will  be  represented 
by  6 x 1 column  matrices 


~ * Nll’  N22’  N12’  Mll*  M22’  M12  ^ 


S = { ell’  e22 


’ 2e 1 2 ’ kll»  k22 ’ 2k12  } 


The  following  elastic  law  is  assumed 
s = E (e  - e") 

where  the  elastic  matrix  E is  the  usual  shell  stiffness 

matrix  relating  the  membrane  forces  N„  and  the  membrane 

strains  e..,  and  the  moments  M..  and  the  curvatures  k.. 

ij  io 

(its  size:  6 x 6). 


The  associated  flow  rule  is  assumed  for  the  plastic 
strain  rates: 


if 


e 


II 


A 


3F 

3s 


F = 1 


and 


(31) 

<‘3s; 


s > 0 


if 


e"  =0 


,3F  T • 

F < 1 or  s < 0. 


(21) 

(22) 


(23) 


■ i 
, * 


a ■ 


The  symbol  3F/3s  stands  for  the  column  matrix 


' . 


and  superscript  T indicates  the  transpose. 


The  parameter  A in  Eq . (24)  can  be  eliminated  with  the 

aid  of  the  condition  F = 1,  or 


F - €)T  2 + (f*)T  2*  = 0 


where 


s * = { 0 , 0,  0,  M n\  M J ) 


II  = { o o 0 

3s*  1 u’  ’ u’  3M*  ’ 3M*  ’ 3M  * 5 


both  being  column  matrices.  Following  a,  by  now,  routine 


procedure,  one  obtains 


,3E\  E e 
3s  ' ~ " 


,3F.T  . 3F  ,3F  ,T  . 3F 

(3J>  E 35  - A 37 


where 


Mn  F 

A - 6(1  - M T2 

0 FH 


Substitution  of  Eq . (24)  into  Eq . (23)  yields 


2 = 5(S  " x H) 


or,  with  A as  in  Eq.  (30), 


s = D e 


where  the  elasto-plastic  tangent  stiffness  D is  given  by 

/!!) T E Al 

D = E[1 — — S y 1 

~ ~ (il)  E iF  _ (II  ) A II 

t'3s;  ~ 3s  ('3s*;  A 3s 


It  scould  be  noted  that  the  stress-strain  relation, 


Eq . (33)  is  of  rate  type.  A proper  numerical  procedure 
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should  be  used  in  evaluating  Eq.  (33)  for  finite  increments 
of  strain,  A e,  and  stress  A s . This  requirement  is,  of 
course,  in  force  for  any  flow  type  theory  of  plasticity. 
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V 


EXAMPLES 


•V 


In  order  to  test  the  present  theory,  the  effect  of  some 


typical  loading  histories  has  been  applied  to  the  following 


shell:  thickness  = 1 in,  Young  modulus  = 29  x 10  psi, 


Poisson  ratio  = 0.3,  yield  stress  in  uniaxial  tension, 


Oq  = 30  x 10  psi. 


Figure  4 shows  moment  - curvature  relation  for  bending 


in  one  plane;  the  curvature  k1 n increases  from  zero  to 

X Ju 


-2  -2 
0.4  x 10  then  decreases  to  -0.4  x 10  and  increases  again 


to  0.4  x 10  . The  results  of  the  present  theory  are  shown 

with  continuous  line.  For  the  same  strain  history,  the 
computations  have  been  performed  with  the  use  of  through-the- 
thickness  integration  (trapezoidal  integration,  21  points 
through  h);  the  corresponding  results  are  shown  with  broken 
line  in  Fig.  4. 


Figure  5 contains  a similar  ca»e  of  bending  in  one  plane 


except  that  the  maximum  and  minimum  values  of  the  curvature 


,-2 


k^^  are  double  of  those  in  Fig.  4,  i.e.  + 0.8  x 10  and 


- 0.8  x 10  , respectively.  It  is  evident  that  further  in- 


crease of  maximum  and  minimum  of  k^  would  not  bring  any  new 


aspects  of  the  moment  - curvature  relation,  since  both  curves 
approach  the  values  M^  = +M?  for  larger  | kn  1 | . 


'll 


The  interaction  of  membrane  forces  and  bending  moments 
is  obviously  important  in  various  problems  of  shell  analysis. 
Figures  6 through  9 show  the  effect  of  interaction  between 


and  M^.  The  loading  histories  in  these  figures  are:  first, 


20 


CX10"1  ) 


Fig.  A Moment  ve.'sus  curv-ture  for  bending  in  one 

plane ; k ^ ^ ^ j ^2  j ” ^12  ~ ^ * ^11  — ^22  ~ ^1*  — ^ * 

continuous  line:  present  theory,  broken  line:  through 

the- thickness  integration 


<X10_i  > 


Fig.  5 Moment  versus  curvature  for  bending  in  one 


plane;  =f  0 


k22 


12 


= 0. 


'11 


= e 


22 


= e 


12 


continuous  line:  present  theory,  broken  line;  through 


the-thickness  integration. 


(XlCf*  ) 


Fig.  6 Moment  versus  curvature  for  bending  and 

extension  in  one  plane;  kn 1 f 0,  e . = 0.25  en, 


<X10  * > 


Fig.  7.  Moment  versus  curvature  for  bending  and 
extension  in  one  plane;  k,,  $ 0,  e, , = 0.50  e„  , 


aio"1  > 


Fig.  8 Moment  versus  curvature  for  bending  and 

extension  in  one  plane;  k ..  ^ 0,  e . = 0.75  en, 


<X10”1  ) 


Fig.  9 Moment  versus  curvature  for  bending  and 

extension  in  one  plane;  k1 1 { 0,  e11  = 1.00  en, 


* V-.  i 


- ■ '-■'  ■-  - 


the  strain  component  e^^  is  increased  to  e^^  = 0.25  e^ 

(Fig.  6),  eix  = 0.50  eQ  (Fig.  7),  e^  = 0.75  eQ  (Fig.  8) 

and  e^  = 1.00  e^  (Fig.  9).  Then,  with  e^  kept  constant, 

the  curvature  is  varied  through  the  cycle  from  0 to 

-2  -2 

+ 0.8  x 10  and  to  - 0.8  x 10  . In  Figs.  6 through  9, 

the  results  of  the  present  theory  are  shown  in  continuous 
line,  with  the  results  of  through-the-thickness  integration 


in  broken  line, 
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VI  CONCLUSIONS 

A theory  of  elas to-plas tic  behavior  of  plates  and  shells 
has  been  presented  in  terms  of  the  membrane  forces  and  moments 
and  the  strains  and  curvatures  of  the  middle  surface.  The 
structure  of  this  theory  is  analogous  to  the  classical  theories 
of  plasticity  of  solids.  It  consists  of  a yield  condition, 
a strain  hardening  rule,  and  a flow  rule.  The  concept  of  the 
yield  surface,  as  known  in  the  classical  plasticity,  exists 
here  in  the  stress  space  of  points  (N^ , N^2»  N.^>  Mll’  ^22’ 

M12*  M21)* 

An  examination  of  the  test  cases  presented  in  this  paper 
indicates  that  the  accuracy  of  the  results  of  the  present 
theory  will  probably  be  acceptable  in  a large  number  of  en- 
gineering applications. 

There  exists  always  the  possibility  of  further  optimization 
of  the  accuracy  of  this  theory.  It  can  be  achieved  by  adjusting 
the  values  of  the  parameters  a and  0,  by  introducing  integer 
or  fractional  powers  of  the  terms  (1  - F^)  and  Fg/F^  in  the 
strain  hardening  rule  (Eq.  16)  etc.  Finally,  some  thought 
should  be  given  to  comparing  the  predictions  of  approximate 
computations  not  to  some  other  theoretical  results  (even  if 
they  are  "exact")  but  to  realistic  experimental  data. 


. ’/  - 
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